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Exciton-polariton propagation in a quantum well, under centre-of-mass quantization, is computed 
by a variational self- consistent microscopic theory. The Wannier exciton envelope functions basis set 
is given by the simple analytical model of ref. [1], based on pure states of the centre-of-mass wave 
vector, free from fitting parameters and "ad hoc" (the so called additional boundary conditions- 
ABCs) assumptions. In the present paper, the former analytical model is implemented in order 
to reproduce the centre-of-mass quantization in a large range of quantum well thicknesses (5aB < 
L < oo). The role of the dynamical transition layer at the well/barrier interfaces is discussed 
at variance of the classical Pekar's dead-layer and ABCs. The Wannier exciton eigenstates are 
computed, and compared with various theoretical models with different degrees of accuracy. Exciton- 
polariton transmission spectra in large quantum wells (L ^ as) are computed and compared with 
experimental results of Schneider et al.^ in high quality GaAs samples. The sound agreement 
between theory and experiment allows to unambiguously assign the exciton-polariton dips of the 
transmission spectrum to the pure states of the Wannier exciton center-of-mass quantization. 

PACS numbers: 78.67.De Quantum wells -71.36.+C Polaritons -71.35.Cc Intrinsic properties of excitons; 
optical absorption spectra 
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I. INTRODUCTION 

Exciton-polariton propagation in mesoscopic super- 
structures has shown new interesting physical phenom- 
ena. Super-radiant effects in multi-quantum wells^— , 
polariton propagation in uniaxial mesoscopic crystals^, 
soliton-polariton and shock waves propagation^, polari- 
tons condensation in DBR cavities^, are present in 
the recent literature. All these phenomena are usu- 
ally observed in exciton-polariton propagation experi- 
ments performed in optically non-local meso-structures. 
Among these, DBR micro-cavities and multi-quantum 
wells under Bragg condition are usually studied as 
paradigmatic systems showing strong radiation-matter 
interactioni^Tii^. For instance, periodic N-quantum wells 
are usually well suited for studying super-radiant exciton- 
polariton propagation, that crucially depends on the op- 
tical properties of theirs 2N-interfaces^'^. Moreover, if 
no-overlapping between the exciton envelope functions, 
localized in adjacent quantum wells, are present, the so 
called dead-layer model was largely used in order to de- 
scribe the optical properties of the 2N-interfaces of the 
systen>i^~— . Recently, a large dead layer effect was also 
observed in colloidal ZnO nano crystals, and this surface 
optical passivation should be responsible on shorten the 
exciton lifetime in this kind of systems 

In his pioneering paper^^, Pekar introduced the con- 
cept of exciton "dead-layer" (or extrinsic dead-layer) in 
order to justify the three-layers model of the electric 
dipole polarization necessary to reproduce the reflectance 
lineshape of Wannier exciton spectra in a semi-infinite 
semiconductor. In a subsequent paper—, Thomas and 
Hopfield discussed the microscopic basis (interface struc- 



tural disorder, chemical disorder, image potential, impu- 
rities etc. . . ) of this classical effect. Looking forward 
a self-consistent microscopic computation of polariton 
propagation in semi-infinite samples in the semiclassical 
framework, D'Andrea and Del Sole^^ introduced an "in- 
trinsic dead layer" (or transition layer) by imposing the 
minimization of the exciton envelope function mismatch 
at the vacuum/semiconductor surface, expanded in an 
hydrogenic basis set. In the same paper an analytical 
approximation of the semi-infinite exciton envelope func- 
tion, that exactly accomplishes the so called no-escape 
boundary conditions (NEBCs) at surface plane (Z=0) of 
the sample, namely: ^/^(^e = 0) = ^^(^^ = 0) = 0, was 
also given for computing the optical response at semicon- 
ductor band edgeS^^: 



^K{r,R)=NK 



(1) 



where f and R are the coordinates of the electron-hole 
relative and center-of-mass motion respectively, (pis(r) is 
the n=l hydrogenic wave function of the relative motion, 
K is the wave vector of the centre-of-mass along Z-axis 
(with Z > 0), is the corresponding in-plane wave vec- 
tor, A = — p~l^ is the exciton reflection amplitude on 

Z=0 vacuum/semiconductor surface (|A|^ = 1) and P is 
an average coefficient of the evanescent waves due to the 
higher energy hydrogenic wave functions (n > 1) taken 
close to the continuum states (n oo)^^?^. Moreover 



the exciton energy is: Ek = — i^*(aB) + ^ 

where R*{aB) is the 3D Rydberg, M and /i are total 
and reduced mass respectively, and is the 3D Bohr 
radius of the exciton. The transition layer 1/P = ^oo 
is strongly dependent on the center-of-mass wave vector 
value K along Z-axis [25], and coincides with the zone. 
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close to the surface plane Z=0, where the exciton enve- 
lope function ^^(^e; ^h] r) ^ 02&i2£. 

Notice that the transition layer of eq.(l) is a quan- 
tum mechanical quantity derived self-consistently by a 
variational principle, and due to the distortion of the ex- 
citon envelope function at surface, where centre-of-mass 
and relative motion of the exciton are entangled, while 
Pekar's dead-layer is a classical quantity characterizing 
the zone of the surface where the electric polarization is 
zero (additional boundary condition- ABC), and due to 
the selvedge layer present in real samples^^. Therefore, 
since Pekar's dead-layer is usually determined as fitting 
parameter, its value can take into account not only the 
effect due to the boundary conditions imposed on the 
exciton envelope function (intrinsic dead-layer), but also 
other physical effects present in real surfaces^^ (extrinsic 
dead- layer). 

The former analytical model of eq.(l) was extended to 
the slab geometry^ where the transition layer and the ef- 
fective Bohr radius were determined by a variational min- 
imization of the first momentum of the exciton Hamil- 
tonian. Also in this case the average level was chosen 
close to the energy of the hydrogenic continuum states 
(n oo)^. Notice, that the former choice, not strictly 
necessary for the consistency of the theory, induces some 
inconsistency as pointed out by the authors of ref. [27,28]. 
In the present paper a different choice, not affected by the 
former problem and that nicely reconciles^^ the transition 
layer in semi-infinite solids and in slabs, is suggested. 

Exciton envelope functions in quantum confined sys- 
tems (wells, wires and dots) are usually computed by 
variational minimization of the exciton energy in a so 
called "ABC free" theoryS^. 

Notice, that in the large quantum well limit {L^ as) 
a very accurate numerical minimization is necessary in 
order to determine a sensible dead-layer value. In fact, 
we remind that the interband transitions at semicon- 
ductor band edge are of the order of some electron- 
volts in the visible range of energies, and the Wannier 
exciton energies are close to the hydrogenic Rydberg: 
R*{aB) = ^-kr (of the order of tens of meV), while 
the centre-of-mass energy of the exciton is given by: 

Ecm{L ^ a^) ~ ^ (x) (of the order of tens of /ieV). 
Finally, the sensitiveness of the optical spectrum to the 
dead-layer effect can be given by the centre-of-mass en- 
ergy difference with and without the dead-layer effect, 
namely: AE{d) = E{L-2d)-E{L) ^ ^ {if ^ (from 
about tens of fieV to zero), where AE{d)/E{L) ^ M/L. 
Therefore, many orders of magnitude^^ in the energy ac- 
curacy should be necessary for computing the Wannier 
exciton dynamics under the centre-of-mass quantization. 

Notice, that the transition layer affects the exciton 
centre-of-mass momentum (spatial dispersion effect), and 
gives additional light waves propagating into the sample, 
therefore its effect on the optical properties cannot be ne- 
glected neither in the presence of a large energy broaden- 
ing. Even if it has a rather negligible effect on the energy 



position of exciton-polariton peaks in transmission spec- 
tra, we expect that it should have a non negligible effect 
on the line shape analysis, and its role should be even 
enhanced in the time resolved optical spectra in multi- 
quantum well systems^i^. 

The aim of the present work is two- folds: i) a non- 
adiabatic Wannier exciton envelope function, obtained 
as a generalization of the analytical model of ref.[l], 
is solved by a variational method, and compared with 
other model calculations at different degree of accuracy, 
namely: a more accurate variational expansion in e-h 
subband products^^i, an exact adiabatic solution^ and 
the heuristic so called "hard sphere" model. Notice, that 
the former extension is necessary in order to study the 
limit of validity of the analytical model itself, and it be- 
comes also mandatory when a complete one exciton basis 
set must be used, as in non-linear optical computation^, 
ii) In order to check the ability of the former analytical 
model in reproducing the exciton centre-of-mass dynam- 
ics, the self-consistent optical transmission spectra in a 
single quantum well is computed and compared with ex- 
perimental results of Schneider et al.^ performed in high 
quality GaAs quantum wells. A sound agreement is ob- 
tained for high quality quantum wells both for intensity 
and phase measurements of the optical transmission spec- 
trum with our analytical model for exciton masses ratio 
from the positronium limit (/i/M = 0.248) to the hydro- 
genic one (/i/M = 0.111); moreover, it allows to unam- 
biguously assign the higher energy peaks, present in the 
spectrum of ref. [2], to the pure states of the centre-of- 
mass quantization. 

In Sect. II the analytical Wannier exciton function 
model of ref.[l] is revisited and generalized. We here 
remind all the approximations involved in order to facili- 
tate the generalization of the analytical model at an order 
higher than the Is hydrogenic wave function as explicitly 
reported in the appendix A. 

In Sect. Ill the optical response of a single quantum 
well is computed and compared with the experimental 
results of Schneider et al^^. 



II. WANNIER EXCITON IN A SLAB: 
VARIATIONAL ANALYTICAL MODEL 
REVISITED AND IMPLEMENTED. 



Let us consider a Wannier exciton of energy E, con- 
fined in a single quantum well of thickness as^ and 
clad between two infinite potential barriers [—L/2 < Z < 
L/2). The exciton function, in a two band model and in 
effective mass approximation, can be expanded in a com- 
plete basis set of hydrogenic eigenstatesi: Hf (pnim {r) = 



3 



*K(r,i?) = 4'£(f,Z) e 



lK/,-R,, 



/Vs 



nlm 



(2a) 



(2b) 



The modulus of the exciton cent re- of- mass w ave 
vector along Z-axis is: Kn = J K'^ — K'j j = 



"2M 



■^^11 where Kji is the wave vector 



along the (x,y) plane. Now, the sum on the hydro- 
genic basis set of eq.(2) can be separated in a finite ( 
for E > Sn ^ h^K'j i/2M) and an infinite sum (with 
E < h^K'j I /2M)^ that, at normal incidence config- 
uration {K 1 1 = 0.0), give the following exciton envelope 
function: 

^E{r.Z)= (3) 

= ^ Wnlm COS {KnZ) + pnlm siu {KnZ)] 
nlm 

nlm 

1/2 

where Pn iKn {Sn - E)] and anlm = Clnlm + 

^nlmi Pnlm — ^ ij^nlm ^nlm^ - 

The analytical approximation of the former wave func- 
tion is obtained by adopting the procedure of ref. [31]. 
First of all, we restrict the finite sum to only one hydro- 
genic function: the Is (n=l) state; the general case will 
be discussed in appendix A and specialized for n=2. 

The computation for a symmetric quantum well pro- 
ceeds by defining the operator: M(z —z; Z —Z) 
that, applied to the exciton envelope function of eq.(3), 
gives the relationship: R ^E(f^Z) = ±^£;(r, Z). From 
this relation two sets of exciton envelope functions are 
generated for even and odd symmetry respectively: 



%{r,Z) = a, cos {K,Z)ifiJr)- 



(4a) 



n>l 



2 ^ [a„im cosh (PnZ) ifil^^ (r) - Pnim sinh (PnZ) (f)] 



*°(f,Z) = I3i Sin {K,Z) if Jr)- 



(4b) 



n>l 



2 ^ [anlm cosh {PnZ) (f^^^ (r) - f3nlm siuh {PnZ) (f^^^ {f)] 
nlm 

where the coefficients ai and /3i are: ai = aioo and 
A = /?ioo. 



In order to perform a smart truncation of the former 
series expansion, let us consider N-hydrogenic functions 
(with N oo), and substitute the infinite sum of vir- 
tual states with an (N-l)-times degenerate state, located 
at an energy value much higher than E. Notice that the 
former approximation can be obtained by substituting 
Pn>i with an average value P = 1/(5, where S is the 
transition layer of the Wannier exciton^ . Therefore the 
former " one transition layer" approximation of Wannier 
exciton envelope function is based on the substitution of 
the infinite sum of the hydrogenic basis set in eqs. (4a) 
and (4b) with four analytical terminators that accom- 
plish the no-escape boundary conditions (NEBCs) at the 
well/barrier interfaces. In fact, by imposing the NEBC: 

^E{Ze = ±L/2,ZhJ) = '^E{Zh = ±L/2,Ze,p) = 0.0 WC 

can determine, as shown in refs. [1,31], the following four 
rational functions: 

n>l 

Fee(r) = 2 ^ {anlm/ai)(pl,^{r) fee{z) (Pls{r) 
nlm 
n>l 

Foo{r) = 2 ^ {anim/ai) (pl,^{r) = foo{z) (pis{r) 

nlm 
n>l 

Feoir) = 2 J2 il3nlm/Pl)^:,J^ = feo{z) Vu{r) 

nlm 
n>l 

F„e(fO = 2 J2 {Pnlm/Hl) <,„(r) = foe{z) </?i,(r) 

nlm 

that give the requested terminators as shown explicitly 
in refs. [1,31] (and also in the appendix A see eqs. (5a)- 
(5)b). This approximation gives rather accurate numeri- 
cal results when a large energy gap is present between the 
travelHng Is state and the virtual higher energy states, 
as will be carefully checked in the calculation and fur- 
ther discussed in this section. With the use of the for- 
mer terminators the exciton envelope functions assume a 
compact analytical form: 

^tnir, Z) = N^icosiKrr^Z) + cosh(P Z) /f™)(^) (5a) 
-sinh(PZ)/f™)(2)]^i,(r) = 
= iV^ g'^iz, Z) ^i,(r) for m = 1, 3, 5, . . . 

<il,{r, Z) = N^[sin{KrnZ) + cosh(PZ) /i™)(^) (5b) 
-sinh(PZ)£™)(^)]^i,(r) = 
= iV° g-^iz, Z) ifuir) for m = 2, 4, 6, . . . 

where Km = Ki{m) with m the centre-of-mass quan- 
tum number for even an odd symmetry, and ^is{r) = 

p-r/c 



/vTra^ is the Is hydrogenic function of the relative 
motion. 

Finally, from the continuity of the exciton wave func- 
tion ^^^\z 0+) = '^^^\z 0") and its first derivative 



dz 



2=0+ 



dz 



for i=e,o at z=0 surface of the 
relative motion, the relationships for the center-of-mass 
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FIG. 1: (Color online) Exciton energy as a function of effec- 
tive Bohr radius (a) and transition layer (S) parameters in a 
AlGaAs/GaAs(001) quantum well of tickness L=250nm. In 
the inset the minimum exciton energy is clearly shown. 



quantization are obtained. Notice, that, due to the sym- 
metry properties of the exciton envelope functions, the 
two former relations are reduced to the condition of ap- 
proaching the z-surface with zero first order derivative. 



namely: 



dF 

e 

dz 



and 



OF 

dz 



0. These con- 



ditions give the centre-of-mass dispersion relations for 
even (m=l,3,5...) and odd (m=2,4,6. . . ) exciton re- 
spectively, namely: 



Kmtg [KmL/2] 
Ptg [KmL/2]- 



^Ptgh [PL/2] = /or m = 1, 3, 5, . . . 

(6a) 

Kmtgh [PL/2] =0/orm = 2,4,6,... 

(6b) 



Obviously, the difference between the former quantiza- 
tion conditions and the simple condition Km = ^ m for 
m=l,2,3,.. is due to the composed nature of the Wan- 
nier exciton and eq. (6a) and (6b) will recover the former 
simple condition in the limit value m ^ oo. 

The variational minimization of the first momen- 
tum of the exciton Hamiltonian is computed for 
the lowest even and odd energy states (m=l,2): 
{^m\Hex \^m) / (^m| |^m) = miu, by taking as varia- 
tional non-linear parameters the effective Bohr radius (a) 
and the inverse of the transition layer (P = 1/S). Even 
and odd variational parameters assume rather the same 
value for large quantum wells {L ^ as) and exactly the 
same value in the bulk limit (L — > oo)r^. 

Since the dynamics of an exciton, perfectly confined 
in a quantum well, strongly depends on the well thick- 
ness (L/a^), let us consider the following different zones 
of quantum well dimensions, namely: i) a zone of very 
large quantum wells (L ^ as) or bulk limit L ^ oo, 
where both the variational parameters assume their bulk 
values, and therefore the exciton energy converges to 
the "hard sphere model", namely: Em = —R*{aB) + 



(L^2d) m=l,2,.... In this zone the tran- 
sition layer assumes its saturation value {5 = 5l^oq)'^^' 
in numerical agreement with that obtained in a semi- 
infinite sample^^'^^. ii) The centre-of-mass quantization 
zone (lOttB ^ ^ ^ 30aB), where the Bohr radius changes 
smoothly and its value remains rather close to the bulk 
one, while the transition layer variation strongly depends 
on the quantum well thickness and dominates the exci- 
ton dynamics from hydrogenic (/i/M=0) to positronium 
behaviour (/i/M= 1/4). Notice, that in this zone the for- 
mer analytival model is well suited for describing the ex- 
citon composite dynamics. Moreover, iii) by decreasing 
the quantum well thickness in the range has ^ L ^ 
lOttB also the Bohr radius value shrinks. In this zone 
strong non-adiabatic effects, due to the entanglement be- 
tween centre-of-mass and relative motion, are present, 
and the analytic model could be implemented by higher 
energy hydrogenic wave function (see appendix A), iv) 
Finally, in the zone where Wannier exciton behaviour 
changes from 3D^2D dynamics, the transition layer and 
the centre-of-mass quantization along Z-axis loose their 
meanings. 

In order to study the Wannier exciton behaviour in 
the centre-of-mass quantization the m=l energy mini- 
mization is computed by adopting the formula reported 
in ref.[31]. In Fig.l the exciton energy, in a rather large 
AlGaAs/GaAs(001) quantum well, is shown as a func- 
tion of the variational parameters value . The physical 
parameter values of the model are the same of ref. [36], 
namely: rUe =0.067mo, rrihh =0.457mo, 715 = 3.71, and 
L = 20aB = 250nm. In this system the heavy hole Bohr 
radius and the Rydberg energy are as =12.4649nm and 
^*(<^B)=-4.1965meV respectively. Therefore, the Wan- 
nier exciton is in the centre-of-mass quantization zone 
with a mass ratio ji/M = 0.111, rather close to the hy- 
drogenic behaviour. 

Notice, that in a scale where the minimum energy with 
respect to the effective Bohr radius is clearly shown, 
the one due to the transition layer is hardly observed 
(see the inset of Fig.l)^. Moreover, the minimized ef- 
fective exciton radius is a=12.4850±0.002nm, and the 
transition layer 5 = 1/P =21.160±0.002nm, while the 
exciton energy for m=l envelope function is £1 =- 
4.17975meV, and the computed centre-of-mass wave vec- 
tor along Z-axis is Km^i = 0.795340 * 10"^a.ii. The 
wave vector and the exciton energy, derived from the 
minimization, are both in rather good agreement with 
those computed by the hard sphere model by taking 
the dead-layer value equal to the minimized transition 
layer {d = S). In fact, the hard sphera model gives: 
Si = -R*{aB) + h^kl/2M = -4.17985meF and Ki = 
7r/{L — 2S) = 0.800525 * 10~^a.u, while exciton wave 
vector computed by zero dead layer value is strongly dif- 
ferent {K = tt/L ^ 0.665012 * 10-^a.u.) from both the 
former values. 

Notice, from eq.(6a) and (6b), that the minimized tran- 
sition layer of the analytical model coincides with the 
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FIG. 2: Effective Bohr radius {a /as) as a function of quan- 
tum well thickness (L/aB) computed for five different masses 
ratios (/x/M :0.050; 0.094; 0.142; 0.219; 0.249). 



dead-layer {S = d) for m=l center-of-mass wave vector, 
while increasing the m- value (m > 1) the dead layer value 
decreases {5 > d)^ till the limit value m ^ oo where the 
centre-of-mass wave vector is Km = mir/L in correspon- 
dence of negligible dead-layer value {d 0) (see also Fig. 
5 and 6). 

In conclusion, the lowest energy exciton state (m=l) 
of our analytical model is in good agreement with the re- 
sults obtained by the heuristic hard sphere model when 
the dead-layer is substituted by the corresponding mini- 
mized transition layer value. Moreover, the sound agree- 
ment with the "exact" adiabatic solution of ref. [29] is 
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FIG. 3: Transition layer values {6 /as) as a function of quan- 
tum well thickness (L/as) computed for five different masses 
ratios {fi/M :0.050; 0.094; 0.142; 0.219; 0.249). 



FIG. 4: Saturation value of the transition layer (Sl^oo/clb) 
as a function of mass ratio /x/M. 



fully discussed in ref. [31], and will not be reported here 
again. 

The transition layer of the non-adiabatic exciton enve- 
lope function of eqs. (5a) and (5b) introduces an entan- 
glement between relative and centre-of-mass motion of 
the exciton, that is the very reason that makes needless 
the ABC imposed to the Maxwell equations, and that 
causes many electromagnetic waves propagation (spa- 
tial dispersion effect) also in semiconductor with cubic 
simmetr y^^i^^ . 

A systematic analysis of the effective Bohr radius and 
of the transition layer as a function of the quantum well 
thickness, and for different values of the mass ratio /i/M 
are shown in Fig. 2 and Fig. 3 respectively. The depen- 
dence of the effective Bohr radius on the well thickness 
(Fig. 2) is a strong function of the mass ratio ji/M. In 
fact, it reaches its saturation limit at about 10 for 
positronium, while the hydrogenic limit converges to the 
correct bulk value for well thicknesses larger than 30 a^, 
and shows a small increment (less than 1/100) with re- 
spect to its bulk value at 7.5 a^. Notice, that also the 
transition layer value reaches its saturation for L?^ 10 
in the positronium, while well thicknesses as large as 
30 are necessary in order to obtain the hydrogenic 
limit (Fig.3). 

The asymptotic values of the transition layer for very 
large quantum wells (L oo), as a function of mass ra- 
tio, are shown in Fig. 4. We observe that the transition 
layer is as large as four-times the Bohr radius in the limit 
of hydrogenic masses ratio {rrih ^ 1836me ^ 5 ^ 4aB), 
while it is a bit greater than half the Bohr radius in 
the positronium limit (M = Aji ^ 5 ^ O.TaB)] its be- 
haviour is also in agreement with the exact adiabatic 
exciton envelope function^^ as shown in a previous com- 
putation for quantum well in the range of thicknesses 
lOa^ ^ L ^ 20aB (see Fig.3 ref. [31] where the two 
curves must be exchanged). We would like to remind 
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FIG. 5: Heavy hole exciton energy as a function of centre-of 
mass quantum number m in a AlGaAs/GaAs(001) quantum 
well of thickness L=250nm. The energy is computed by full 
theory (dots), analytical model (crosses), hard sphere model 
with (solid curve) and without dead layer (dashed curve). 

that a former evaluatior>2^, derived from the energy bal- 
ance equation: -i?*(a) + jM-^m ~Jm (i) ' 
continuum limit of the hydrogenic set of states gives not 
a correct behaviour of the transition layer— as a function 
of exciton masses ratio (see also ref.[31] and the discus- 
sion therein). 

Notice that exciton transition layer of Fig. 4, computed 
in the limit of slab thickness L ^ oo, should be equal to 
those in semi-infinite sample of eq.(l). Therefore, the 
former choice can reconcile the transition layer effect in 
semi-infinite samples^^^^^ and in slabs^, without intro- 
ducing a further "ad hoc" dead layer (added to the tran- 
sition layer) as hypothesized in ref.[34]. 

We have observed before, that in the analytical model 
of eqs. (5a) and (5b) the most intriguing approxima- 
tion is the hydrogenic series truncation adopted. In the 
present paper, the correct energy location of the average 
virtual state^^ is obtained by comparing the analytical 
model with the full theory of ref. [32]. 

In Fig. 5 and 6 the heavy and light hole exciton ener- 
gies, computed by the full theory of ref. [32] and by the 
former analytical model for a quantum well of thickness 
L=20aB=250nm are given as a function of an energy or- 
der parameter m that in the case of our analytical model 
coincides with the centre-of-mass quantum number. The 
physical parameters adopted in the calculations are the 
same of ref. [36], and the minimized parameters are: 
aih = 19.977 ±0.002nm and ahh = 12.483 ± 0.002nm for 
light and heavy hole effective Bohr radii respectively, and 
6ih = 14.277 ± 0.002nm, Shh = 21.133 ± 0.002nm for the 
corresponding transition layers. 

We observe, from Fig. 5 and Fig. 6, that the contribu- 
tion of n=2 hydrogenic function becomes important for 



heavy hole exciton states with the centre-of-mass quan- 
tum number greater than m>12 corresponding to the 
exciton energy ~ 2.3meV greater than the heavy hole 
m=l state, and analogously for light hole states m>8 
(^ S.SmeV). Therefore, the analytical model, also in 
its lowest order of formulation^, should be able to re- 
produce the optical response in this range of quantum 
well dimensions since both the former ranges of energy 
are greater than those shown in the experimental spectra 
(see Fig. 14b) of Schneider et alS^. 

Finally, in Fig. 5 and Fig. 6 also the "hard sphere" 
model energies, computed with and without {d 0) the 
dead layer effect, are shown. The rather sound agree- 
ment between the full theory and the hard sphere model 
with dead-layer effect, for the lower m- values is essen- 
tially due to the values adopted in the present calcula- 
tion, that coincide with the transition layer values ob- 
tained by a variational minimization (namely: S^h = 
1/Phh ^21.1nm and Sih = '^/Pih ^14.4nm) while those 
computed with the hard sphere interpolation equation: 
d = {1 — 2/i/M) as are strongly different from the for- 
mer ones {dhh ^ lOnm and dih ~6.5nm for heavy and 
light hole exciton respectively). 

Now, let us compute the lowest energy envelope func- 
tion of an exciton in the hydrogenic limit {ruh oo) for 
the heavy hole trapped at the site = of a rather 
large quantum well {L = 20aB)- The analytical envelope 
function is: ^m(^;Zo = 0) = Nmgm{z;Zo = 0)(pis{r), 
and, for transition layer limit S ^ {P ^ oo), the 
exciton energy converges to the bulk value. Indeed, in 
this limit, from the eqs. (5a) and (5b) , we obtain 
tg{KmL/2) — oo and tg{KmL/2) = respectively, 
therefore the wave vector quantization of the electron is: 
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FIG. 6: Light hole exciton energy as a function of centre-of 
mass quantum number m in a AlGaAs/GaAs(001) quantum 
well of thickness L=250nm. The energy is computed by full 
theory (dots), analytical model (crosses), hard sphere model 
with (solid curve) and without dead layer (dashed curve). 
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function shape. Therefore for very large quantum wells 
{L > 20a b) the former analytical model should be able to 
describe the optical experimental mesurements of ref. [2] . 

In order to generalize the former model, also in the 
zone very close to the quantum well thicknesses Sa^ ^ 
L ^ lOa^, we will have to extend the finite sum of eq.(3) 
to higher energies of the hydrogenic wave function, as 
discussed in the appendix A. Moreover, since the real 
symmetric matrix Hamiltonian is rather diagonal only 

for large wells, in this range of thicknesses, we must solve 

— > -(-^ — > 

a generalized eigenvalue problem, namely: A = E B 

where Hex |^m) = (^m| |^n) = ^ and the ma- 
trix of the normalized exciton states is quasi-unit real 

symmetric matrix (^^1 |^m) = (^m| |^n) = B. Notice, 
that the generalized analytical envelope function is not 
a pure state of the center-of-mass wave vector, even if it 
still remains under the one-transition layer approxima- 
tion. 

In conclusion, we have pointed out that the analytical 
model is essentially based on a smart truncation of the 
hydrogenic basis set expansion, and it must be imple- 
mented (see appendix A) to correctly take into account 
the Wannier exciton dynamics in all the range of thick- 
ness till the semi-infinite limit (Sa^ ^ ^ ^ oo). 
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FIG. 8: Trapped exciton confinement function 
\gm=i(ze; Zh = 0)\ of AlGaAs/GaAs(001) quantum well 
of thickness L=250nm as a function of electron Z-coordinate 
with transition layer S = 21.16nm. 



Km = ^ m for m=l,2,3,... . In appendix B the first 
momentum of the lowest energy state of the trapped ex- 
citon Hamiltonian in a large quantum well is given in 
an analytical form. The bulk exciton energy can be ob- 
tained also by a variational minimization as a function 
of effective Bohr radius a = as = 10.8711 ± 0.002nm 
and the effective Rydberg is Ryd=-4.81173meV. In Fig. 7, 
the confinement function gi{z] Zq = 0) is computed as a 
function of the electron-hole relative motion along z-axis. 
Notice, that, at variance of a simple cosine function that 
describes electron confinement in a 2D quantum well2i, 
this function shows a rather flat shape, except at the 
interfaces (Z ^ ±L/2), and, as fully discussed in ref. 
[27,28], this behaviour is a fingerprint that the present 
analytical model is well suited for describing the centre- 
of-mass quantization. In Fig. 8 we compute the former 
system by imposing a transition layer value different from 
zero (b =21.133nm); also in this case the shape remains 
rather smooth, but the transition layer effect is clearly 
shown at the boundaries of the well. 

In conclusion we would like to remind that the opti- 
cal response of exciton-polaritons is strongly dependent, 
not only on the exciton energy, but also on its envelope 



III. OPTICAL RESPONSE OF A SINGLE 
QUANTUM WELL UNDER THE 
CENTER-OF-MASS QUANTIZATION. 

In this section, the exciton-polariton propagation in 
a single quantum well under centre-of-mass quantiza- 
tion is computed by the former analytical exciton model 
in a large range of quantum well thicknesses (Sa^ < 
L < 20aB ), and compared with the full theories of 
refs. [2,31,36], and with the heuristic hard sphere model. 
Moreover, we will compute self-consistently the polari- 
ton propagation in a slab in order to reproduce, without 
introducing fitting parameters (except the homogeneous 
non-radiative broadening), the experimental transmis- 
sion spectra in the high quality GaAs samples of Schnei- 
der et al.^ 

First of all, let us consider a GaAs quantum well of 
thickness L = has = 62.5nm, that is in the transition 
zone, as discussed in the former section. The parame- 
ter values, adopted for the calculation, are taken from 
ref. [36], namely: rrie = 0.067mo, rrthh = 0.457mo, back- 
ground refraction index = 3.71, energy gap Egap = 
lA2eV, and non-radiative broadening r=0.03meV. The 
electric dipole moment is replaced by the Kane's energy^ 
{Ek =22.71eV for GaAs), that is conceptually more suit- 
able than the electric dipole moment for describing the 
radiation-matter interaction in translational periodic sys- 
tems. 

The transmittance spectrum is computed by a self- 
consistent solution of Schroedinger-Maxwell equations in 
the effective mass approximation, as reported in ref. [31], 
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FIG. 9: Optical absorbance (1-T) of the heavy hole exci- 
ton as a function of photon energy {{E — Egap) / Ryd) in Al- 
GaAs/GaAs(001) quantum well of thickness L=62.5nm. 



where all the analytical formula necessary for the cal- 
culation are given explicitly. The Fabry-Perot effect is 
suppressed in the spectrai^s^ by using the same back- 
ground dielectric function in whole the space; this allows 
to observe exciton polariton peaks free from the interfer- 
ences due to the Fabry-Perot oscillations^. In Fig. 9 the 
absorption spectrum is shown for the lowest five polari- 
ton states; the energy positions of the optical peaks are 
labelled with the same quantum number m of the cor- 
responding centre-of-mass exciton state (m=1^5). For 
this quantum well thickness (L=A/4) the even and odd 
exciton-polariton states show almost comparable intensi- 
ties. 

The phase of the transmittance amplitude is shown in 
Fig. 10. Notice that the maxima of the absorption are 
in correspondence of the inflection points of the phase 
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FIG. 10: Phase of the optical absorbance of the heavy hole 
exciton as a function of photon energy [[E — Egap) / Ryd) in 
AlGaAs/GaAs(001) quantum well of thickness L=62.5nm. 




FIG. 11: Heavy hole exciton energy as a function of centre- 
of-mass quantum number m in AlGaAs/GaAs(001) quantum 
well of thickness L=62.5nm. The energy is computed by full 
theory (dots), analytical model (crosses), hard sphere model 
with (solid curve) and without dead layer (dashed curve). 



curve; therefore, there is a complete correspondence in 
the polariton spectra between phase and intensity. To 
properly reproduce both these quantities is a very severe 
check for the exciton model adopted in the calculation as 
fully discussed in ref. [2] and underlined in the present 
theory versus the experimental results. 

Now, let us compare the exciton energies computed 
by the analytical model with those derived by the full 
theory of ref. [32]. In Fig. 11 these exciton energies 
are shown as a function of the quantum number m 
of the centre-of-mass motion. The minimized effective 
Bohr radius and transition layer of the analytical model 
are: a=12.464±0.002nm and (5=1/P=10.228±0.002nm 
respectively. The exciton energy computed with the sim- 
ple hard sphere model are also reported for vanishing 
dead layer and for a dead layer as large as the minimized 
transition layer value. We note a rather large energy 
difference (about l,23meV for m=4) between analytical 
model and full theory in correspondence with higher en- 
ergy states (m>4), that rapidly decreases for the lower 
energy states; for instance for m=3 it drops to O.OmeV 
and for m=2 is about 0.2meV. The last value is rather 
acceptable if we take into account that the non-radiative 
broadening is about Tnji = 0.5^1.0 meV in tipycal sam- 
ples and Tnh ^ 0.25meV in high quality quantum wells. 
Obviously, the large discrepancy observed for the high en- 
ergy states is not surprising if we take into account that 
the exciton envelope functions of the analytical model, 
labelled with the quantum number m, are pure states of 
the centre-of-mass, while the exciton states of the full 
theory are a general superposition of states with differ- 
ent centre-of-mass wave vectors and are expanded on a 
rather complete set of the relative motion wave functions. 

In conclusion, in this thicknesses range of values the 
lowest exciton states are in sound agreement with those 
obtained by the full theory, while for describing higher 
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FIG. 12: Heavy hole exciton energy as a function of centre- 
of mass quantum number m in AlGaAs/GaAs(001) quantum 
well of thickness L=125nm. The energy is computed by full 
theory (dots), analytical model (crosses) , hard sphere model 
with (solid curve) and without dead layer (dashed curve). 



energy states of the centre-of-mass (m>4) we will have 
to implement the simple analytical model as discussed at 
the end of the former section. Moreover, we observe that 
the full theory gives exciton states located between the 
two curves of the hard sphere model with and without 
the dead layer effect and this leads to the well known pos- 
sibility of recovering the rather correct exciton energies 
by using the Pekar's dead-layer as a fitting parameter. 
Now, let us consider GaAs quantum wells in the zone of 
the centre-of-mass quantization lOa^ ^ L < 20a b- In 
Fig. 12 the exciton energies of the analytical model and 




(E-Egap)/Ryd 



FIG. 13: Optical absorbance (1-T) of the heavy hole exci- 
ton as a function of photon energy {{E — Egap) / Ryd) in Al- 
GaAs/GaAs(001) quantum well of thickness L=125nm: full 
theory (dashed curve), analytical model (solid curve). 



of the full theorys^ in a single quantum well of thickness 
L=10aB=125nm are shown as a function of the centre- 
of-mass quantum number m. The physical parameters 
values are the same of the former sample, and the mini- 
mized variational parameters are: a = 12.560 ± 0.002nm 
and (5 = 1/P = 16.400 ± 0.002nm. We observe that, 
while the effective Bohr radius is very close to the 3D 
Bohr radius, the transition layer value is rather different 
from its saturation value, underlining the strong influence 
of the transition layer effect on the exciton dynamics in 
this range of thicknesses. The calculation confirms a ne- 
glegible difference between the two models for the lower 
energy states (m ^ 4), and an increase of this differ- 
ence till 0.75meV for the m=8 state where higher energy 
(n > 1) hydrogenic states cannot be neglected. In fact, 
exciton states with m>8 drop into the second parabolic 
curve due to the n=2 hydrogenic states, therefore the 
agreement between the analytical model and the full the- 
ory can be improved along the line discussed in the ap- 
pendix A. In Fig. 12 also the exciton energies computed 
by Pekar's model with and without dead-layer are shown 
for sake of comparison. 

In Fig. 13 the optical spectra, computed self- 
consistently with the analytical exciton envelope function 
and with the full exciton model, are shown; due to the 
photon wavelength that is about A ~ 2L the odd exci- 
ton functions (m=2,4,6,. . . ) give a stronger contribution 
than the even (m=l,3,7,. . . ) ones. In fact, in the com- 
puted spectrum of the full exciton theory^^ the intensity 
of m=7 exciton-polariton peak is vanishingly small, while 
it is shifted towards the higher energies for the analytical 
model. We would like to underline that the transmission 
spectrum computed by our full theory is in very good 
agreement, both in energy positions and in line-shapes, 
also with those computed by the full theories of refs. [2] 
and [36], not reported here. 

In order to go a bit deeper in the optical spectrum 
analysis, notice that the double peak, due to the low- 
est exciton energies (m=l,2), is very sensitive to the 
radiation- matter interaction. Moreover, due to the small 
non-radiative broadening value adopted in the calcula- 
tion, we can appreciate the difference in energy between 
the analytical model and the full theory for all the range 
of energies of the polariton spectrum even if theirs val- 
ues remain in any case rather small for m < 8 exciton- 
polariton states. In fact, let us consider the two peaks 
for m=6: they seem very well separated in the spectrum 
but the energy difference (<0.4meV) is lower than the 
optical broadening (see also Fig. 12). 

Finally, we would like to underline that in the range of 
energies reported in Fig. 13 our microscopic model gives a 
more accurate result than that based on Pekar's ABC. In 
fact, in ref. [36] (see fig. 4) the energy difference between 
full theory and Pekar's model for m=6 exciton-polariton 
state is about 0.7 meV, that is two-times our difference, 
and it is greater than the broadening value. Therefore, 
we can conclude that our analytical model can reproduce 
the full theory till m < 6 exciton-polariton state, and 
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that, in any case, it gives a better agreement with the 
full theories than that obtained by Pekar's mode^. 

Now, let us try to reproduce the experimental trans- 
mission spectrum of an high quality GaAs quantum well 
reported in the work of Schneider and al.^ by our analyti- 
cal self-consistent theory. In fig. 7b of ref. [2], reproduced 
here in Fig. 14b, many peaks due to the centre-of-mass 
quantization of both heavy and light hole excitons in a 
large quantum well of GaAs {L = 20aB = 250nm) are 
clearly shown by arrows. We would like to underline that 
we have choosen to reproduce the experimental results of 
Schneider et al.^ because they are well suited for check- 
ing the present analytical model. In fact, the Wannier 
exciton is perfectly confined in a self-sustained quantum 
well, and the Fabry-Perot effect, due to the modulation of 
the background dielectric function, is suppressed by the 
presence of antirefiection coating on both the sample sur- 
faces. Moreover, in a large quantum well the anisotropic 
electron-hole masses induce a very small energy shift on 
the exciton energies, therefore, we can use isotropic exci- 
ton masses in the present range of thicknesses. 

In the former section, we have shown that the analyt- 
ical exciton model is in perfect agreement with the full 
theory^S in a sufficient large range of energies (9meV), 
and this is well suited for reproducing the experimental 
spectrum of ref. [2]. The parameter values adopted for 
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FIG. 14: (a) Optical absorbance (1-T) computed for light 
and heavy hole exciton as a function of photon energy in Al- 
GaAs/GaAs(001) quantum well of thickness L=250nm. (b) 
Experimental optical absorbance of light and heavy hole ex- 
citon of ref. [2] as a function of photon energy in Al- 
GaAs/GaAs(001) quantum well of thickness L=250nm. 



= O.OSOmo, E'Jlp = 1.5105ey, 



the calculation are^: L 
rrihh 0.457mo, ruih 

Ub = 3.71, E^^p = 1.5130eV and Kane's energy Ek = 
22.1\eV\ the minimized parameter values, given in the 
former section (see also Fig. 5 and Fig. 6), underline the 
different behaviours, from positronium for light hole to 
hydrogenic atom for heavy hole, of Wannier excitons. 
The energy shift between heavy and light hole energy 
gaps, due to the residual strain of the sample, is discussed 
and evaluated in ref. [2]. 

The high quality {J^nr ^ 0.25mey) of the samples 
should allow to study the effect induced by the wave vec- 
tor quantization of the centre-of-mass motion at T=2K. 
Indeed, the polariton effect is important only in the case 
where the exciton broadening value is less than a critical 

valuei^, namely: = [^^ig^^j , that for the pa- 
rameters chosen is Fc ^ 0.21meV. Notice that the non- 
radiative broadening value chosen in the present calcula- 
tion, in order to reproduce the experimental spectra line- 
shapes {Tnr = O.lSmeV), is in rather close agreement 
with that given in ref. [2], and, moreover, since F < Fc, 
exciton-photon coupling becomes dominant with respect 
to the exciton acoustic phonon interaction. 

Since the exciton is computed in a simple two bands 
model, the optical response of heavy and light hole 
exciton-polaritons will be shown in the pictures by dif- 
ferent curves but in the same energy scale. 

The absorbtion of light and heavy hole excitons are 
shown in Fig. 14a. We would like to underline that the 
agreement between theory and experiment is very good 
both in the line shapes and in energy positions, as shown 
by the arrows in correspondence of the maxima of the ex- 
perimental spectrum of Schneider^ reported in Fig. 14b. 
In fact, the characteristic double peak of the light hole 
exciton-polariton and the asymmetric line shape of the 
main peak of the heavy hole, due to the interplay of 
polaritonic effect and centre-of-mass quantization of the 
Wannier exciton, are clearly reproduced. This agreement 
is particularly meaningful when we take into account that 
it is obtained for the same parameter values in the pres- 
ence of rather different exciton dynamics. In fact, light 
hole exciton is close to positronium limit (/i/M=0.248), 
while the heavy hole is rather close to the hydrogenic be- 
haviour (/i/M=0.111). Moreover, the analytical model 
allows to unambiguously assign the center-of-mass char- 
acter for any features present in the experimental spec- 
trum. For instance, the peak at m=7, not assigned in 
ref. [2], has even parity character, due to the photon wave- 
length {L ^ A) that gives odd exciton-polariton peaks 
neglegible small in the spectrum. 

Finally, notice that for the heavy-hole exciton polari- 
ton in ref. [2] (see Fig. 10c) the m=ll peak is shifted to- 
ward the lower energy side, with respect to the experi- 
mental spectrum, while the m=13 one is out of the energy 
scale. 

In Fig. 15a and 15b the computed phase of exciton- 
polaritons absorbtion spectrum and the corresponding 
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FIG. 15: (a) Phase of the optical absorption computed for 
Ught and heavy hole exciton as a function of photon energy 
in AlGaAs/GaAs(001) quantum well of thickness L=250nm. 
(b) Measured phase of the optical absorption^ of light and 
heavy hole exciton as a function of photon energy in Al- 
GaAs/GaAs(001) quantum well of thickness L=250nm. 



experimental results of ref.[2] are shown respectively. 
Also in this case the agreement between theory and ex- 
periment is sound both for energy positions and line- 
shapes. Obviously, the long tail shown by the heavy- hole 
exciton-polariton phase towards the lower photon ener- 
gies can not be matched with the light hole phase by 
the simple two-bands model. Moreover, all the features 
observed in the absorption spectrum have a counterpart 
in the phase spectrum, and also the characteristic two 
peaks feature of the light hole exciton-polariton and the 
centre-of-mass peaks at higher energies are reproduced. 

In conclusion, the simultaneous measurements of phase 
and transmission spectra allow to assign the exciton- 
polariton peaks of the experimental spectrum to the pure 
centre-of-mass exciton-polariton quantization. 

We would like to underline that in ref. [2] the au- 
thors compare the optical properties, computed by their 
accurate microscopic exciton model, with a number of 
heuristic models, based on different additional boundary 
conditions (ABCs), in order to conclude that "A simulta- 
neous description of the measured amplitude and phase 
of the transmitted electric field is only possible with a 



full model". The sound agreement between our analyti- 
cal microscopic model, based on the one-transition layer 
approximation, and the experimental results of ref. [2] 
strongly suggest the replacement of the adjective full 
with microscopic. 



IV. CONCLUSIONS 

In the present work the analytical microscopic exciton 
model in a large {L^ as) single quantum well, proposed 
in ref. [1] is revised and implemented in order to repro- 
duce the Wannier exciton centre-of-mass quantization in 
the range of quantum well thicknesses bas ^ L ^ oo. 

The microscopic transition layer effect, observed at the 
well/barrier interfaces in a single quantum well, is dis- 
cussed at variance of the classical dead-layer introduced 
by Pekar^^, and compared with various full theoretical 
models^^^. 

Moreover, this new formulation has allowed to recon- 
cile the Wannier exciton dynamics in large quantum wells 
{L ^ as) and in semi-infinite semiconductor samples 
(L ^ oo) described by the same analytical microscopic 
modeli^. 

Finally, exciton-polariton transmission spectra in large 
quantum wells (L = 20aB) are computed, and compared 
with experimental results in high quality GaAs samples 
obtained by Schneider et al.^. The sound agreement be- 
tween theory and experiment has allowed to unambigu- 
ously assign the transmission spectra peaks to the pure 
states of the exciton centre-of-mass quantization. 
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Appendix A: A generalized variational Wannier 
exciton envelope function 

We generalize the Wannier exciton envelope function 
expanding eq.(2) of the text by taking into account also 
the 2s hydrogenic function. Since 2s is symmetric with 
respect to the z-axis we can adopt the same procedure of 
Section II based on the following three steps. 

1) From the relationship R^£;(r, Z) = ±^£;(r, Z) we 



12 



obtain the even and odd envelope functions, namely: 

n>2 

2 [anim cosh {PnZ) (r) - f3nim siuh {P^Z) ^^^^ (r)] 



and analogously for the odd ones: 

sin(KZ2) sinh(PZi) - sin(i^Zi) sinh(PZ2) 



nlm 



(Ala) 



foo{K) 



feo{K) 



sinh [P{Z2-Zi)] 

(A6a) 

cos{KZ2) cosh(PZi) - cos{KZi) cosh(PZ2) 



sinh [P(Z2-Zi)] 



(A6b) 



n>2 

2 ^ [a„;™ cosh (P„Z) (f) - sinh (P„Z) ^l^^ (f)] 



nlm 



(Alb) 



where ai = aioo, /^i = /^loo; the functions ^^^(f^Z) for 
i=e,o are: 

$^(r,Z) = cos(i^iZ) ^is{r) +0^200 cos(i^2^) V^2s(r) 

(A2a) 

(r, Z) = sin {KiZ) (pi,(r) + ^200 sin {K2Z) i^2s{r) 

(A2b) 

where 0^200 = <^20o/<^i and ^^200 = hoo/ Pi- Indeed, the 
functions of eqs.(la.A) and (lb. A) will be reduced to 
eqs.(4a.A) and (4b. A) respectively for a200 = hoo = 0. 

By adopting the one transition layer approximation an 
average value for the inverse of the transition layer depth, 
Pn>2 = P = 1/^^ is taken and the envelope functions 
assume the following analytical form: 

n(r,Z) = 

= Ne [^|(r, Z) + cosh(PZ)Fee(r) - sinh (PZ) Peo(r)] 

(A3a) 

^Ur.Z) = 

= No [^|;(r, Z) + cosh(PZ)P,,(r) - sinh (PZ) P,e(r)] 

(A3b) 

where = ai and Nq = Pi. 

2) We impose the no-escape boundary conditions 
(NEBCs):^ ^E{ze = ±L/2,Zh,p) = ^ E{zh = 
±L/2,Ze,p) = 0.0 in order to compute the F-functions. 
For the even solution ^|;(r, Z), eq.(A3a), we obtain: 



Fee{r) = fee{Ki) (pis{r) + Q^2 /ee(i^2) ^2s{r) 
Feo{r) = feo{Ki) (pis{r) + Oi2 feo{K2) ^2s{r) 



(A4a) 
(A4b) 



for < z < P/2, Zi{z) = aeZ-L/2 and Z2{z) = -anz- 
L/2. The f-functions for even exciton functions are: 



fee{K) 
foe{K) 



cos{KZ2) sinh(PZi) - cos(i^Zi) sinh(PZ2) 
sinh [P(Z2-Zi)] 

(A5a) 

sin(i^Z2) cosh(PZi) - sin(i^^Zi) cosh(PZ2) 



sinh [P(Z2-Zi) 



3) From the continuity of the exciton wave function 
and its first derivative at the z=0 surface we obtain the 
conditions for the centre-of-mass quantization. Also in 
this case, due to the symmetry properties of the different 
components of the envelope functions, the centre-of-mass 
quantization must accomplish the two following relation- 



ships: 



dz 



Oand ^ 



for the even exci- 



ton function (analogously for the odd one). It is simple 
to check that the centre-of-mass wave vector K follows 
the even dispersion: K^tg[K^L/2] + Ptgh[PL/2] = 
for m=l,3,5,..., and the odd one: Ptg[K^L/2] - 
K^tgh [PL/2] =0 for m=2,4,6,.., according to the even 
{fee{z)Jeo{z)) and odd {foo{z)Joe{z)) f-fuuctious respec- 
tively. The generalized even (m=l,3,5,...) exciton func- 
tion of eq.(A3a), assumes the analytical form: 



^mi,m2 {cOs{Km^Z)(pis{f)^ ^2 COs{Krr 



(A5b) 



(A7) 

^m2Z)(f2s{f) + 
+ COsh(PZ) [fee{KmJ(pis{r) + Q^2/ee (i^ms ) V^25 (0] + 

- sinh(PZ) [feo{KmJ(pis{r) + a2feo{Km2)^2s{r)]} 

and analogously for the odd Wannier exciton function. 

Notice that the inclusion in the analytical model of 
the 2s hydrogenic function not only increses the number 
of variational parameters, but also gives off-diagonal en- 
velope function components in the centre-of-mass wave 
vector. 

The inclusion in the model of 2p hydrogenic functions 
is a bit more tricky, since the finite sum is composed 
of two even components {2px^ '^Py)^ but also of an odd 
one {2pz). Therefore, the generalized envelope functions 
should mix even and odd centre-of-mass wave vectors. In 
fact, let us choose (^,/,C) as the direction cosines of the 
relative motion coordinate with respect to the Cartesian 
axis: (p2p{r) = ^^2p^{r) ^ f ^2py{r) + C^2p,(r). In this 
case eqs. (A2a) and (A2b) become: 

<l>^(r, Z) = cos(L:iooZ) (pis{r) +0:200 cos(L:200^) ^2s{r)- 
+0211 cos(L:21iZ) [i(p2p^{r) + f^2py{r)] + 
+C/^2io sin {K210Z) (p2p, (r) 

(A8a) 

<l>|;(r, Z) = sin {KiqqZ) ^uir) + p2QQ sin {K200Z) (f2s{r)- 
+^211 sin {K211 Z) [£^2p. (r) + f^2py (r)] + 
+Cck^2io COS {K210Z) (p2p, (r) 

(A8b) 

where 0200 = <^20o/<^i,^Ok;2ii = 0211/0^1, P210 = /^2io/<^i 
and P200 = /^20o//^i, P211 = /52ii//^i, <5^2io = c^2io/Pi- 
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Now, let us impose the no-escape boundary condi- 
tions (NEBCs): ^E{ze = ±L/2,Zh,p) = ^E{zh = 
±L/2^ p) = 0.0 in order to compute the F-functions. 
For the even solution ^|;(r, Z) of eq.(A3a) we obtain: 

+Q^21l/ee(i^21l) I ^2pA^) ^ f ^2py{r) + 

(A9a) 

Feo{r) = feo{Kioo) ^ls{r) + Oi2 feo{K200) ^2s{r)^ 
+/^21l/oe(^21l) i(P2pA^) ^ f^2py{r) + 
+Cc^210/oe(i^210)v^2p,(r) 

(A9b) 

Finally, the generalized even (m=l,3,5,...) exciton func- 
tion derivation proceed as before. Notice, that the for- 
mer generalization leads to increase the number of vari- 
ational parameters; in fact, while for n=l we have two 
non-linear parameters (a, (5), for n=2 we have two non- 
linear (a, (5) plus three linear (0^200, <5;2ii, $211) parame- 
ters of minimization. Analogously for odd exciton state 
(m=2,4,6,...), and also in this case the solution embodies 
odd and even centre-of-mass wave vectors. 

Finally, by adopting the former procedure, we obtain a 
set of normalized independent exciton states, and after- 
wards, it is necessary to solve the generalized problem in 
order to obtain an orthogonal set of states as discussed 
in Section II. 



Appendix B: Trapped exciton in a quantum well 

Let us consider a Wannier exciton trapped as a neutral 
atom at the site Rq = (0,0,Zo), whose electron mass is 
me, and the hole mass is taken in the limit rrihh ^ 00. In 
this case the Z-component of the centre-of-mass coordi- 
nate is: Z = Zh = while the coordinate of the relative 
motion is: z = Ze — Zh = Ze — Zq] therefore, the motion 
of the electron along Z-axis is confined in the segment: 
—Zo — L/2 < z < — Zo + L/2. Since the envelope func- 
tion of the trapped exciton has not defined parity, except 
for the site = 0, a general variational trial function 
can be obtained as a linear superposition of even and 
odd exciton envelope functions in order to remove the 
restricted symmetry properties. Here, we adopt a more 
intuitive procedure observing that even and odd exciton 
envelope functions, defined in Sect. II, lost their symme- 
try property if projected into a non symmetric domain 
as the segment —Zq — L/2 < z < —Zq + L/2. In this 
case, m=l envelope function, projected into the segment 
—Zq — L/2 < z < —Zo + L/2, can be used for a vari- 
ational determination of the effective Bohr radius that 
should converges to the 3D Bohr radius for large quan- 
tum well (L >> a^), while the transition layer assumes 
the limit (5^0. Obviously, many physical properties of 



exciton analytical model have to be redefined for trapped 
exciton. In fact, the projected envelope function is: 

^fm{r]Zo) = N,ngm{z]Zo)^is{r) (Bl) 

For m odd (m=l,3,5,) we have: 

gm{z] Zo) = cos(K^Zo) + 

+ COsh(PZ,)/ee(^) - sinh(PZ,)/eo(^) (B2) 

and for m even (m=2,4,6,..) : 

gm{z', Zo) = sin(i^^Zo) + 

+ cosh(PZ,)/,,(^) - sinh(PZ,)/,e(^) (B3) 

where the f-functions are obtained by imposing the 
NEBCs at the surfaces of the weh (Z = ±P/2). For 
< z < L/2 the confinement boundaries are Zi = 
-Zo - L/2 = -L/2 and Z2 = -Zo + P/2 = -z + L/2: 

fee{z) = [cos(L:Zi)sinh(PZ2) - cos(L:Z2) sinh(PZi)] 
feo{z) = [cos{KZi) cosh(PZ2) - cos(L:Z2) cosh(PZi)] 

(B4) 

and for —L/2 < z < the confinement boundaries are 
Zi = L/2 and Z2 = -z - L/2: 

foe{z) = [sin(L:Zi) cosh(PZ2) - sin(i^Z2) cosh(PZi)] 
foo{z) = [sm{KZi) sinh(PZ2) - sin(L:Z2) sinh(PZi)] 

(B5) 

and fij{z) = fij{z)/ sinh[P{Z2 - Zi)] for i , j=e ,0 . 
Remember that the continuity of the exciton envelope 
function and its first derivative at the surface z=0 are 
the same of those given in eq.(5a) and (5b), but in the 
trapped exciton appears the electron quantization con- 
dition at variance of the former case that describes the 
center-of mass motion of Wannier exciton. 

Taking into account that we have two different analyt- 
ical equations for 2: > and 2; < 0, the normalization 
integration is: 



(^mll^m) 



(B6) 



27r 



00 ^ U Z2 \ 

j pdplj dz + j dz \gl{z;Zo)^ls {V?T^ 



and the first momentum of trapped exciton Hamiltonian: 



HI 



-Vl- 



is: 



2/i 

\^m)=- — \{^m\\^m) 



(B7) 



CO 








"1 1 " 


j pdp< 




-P 










[me 


a as 



9li{z;Zo) 



h^gm{z]Zo) 
2me 



(1^1 + 2)^ 



2^9grj 



a dz 
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The n-th derivative with respect to the relative motion 
along z-axis of the confinement functions gm{z] Zq) are: 



- cos(PZo) 



• sin(PZ^) 



(B8) 



dz^ dz^ dz^ 

Taking into account that dZ2 = —dz^ the derivative are: 



Finally, for very large quantum wells the variational en- 
ergy of trapped exciton recover the 3D Rydberg value, 
while the confinement function \gm{z; Zo) \ is an asym- 
metric function except for the hole coordinate Zo = Zh = 
and —L/2 < Ze < L/2 as discussed in the text (Fig. 7 
and 8). 



9fij 
dz 



and 



P/,,cosh [P(Z2-Zi) 



9fij 
dz 



dz^ ~ ^ 



dz^ 



2P 



^cosh[P(Z2-Zi) 



|/sinh [P(Z2-Zi)] 
(B9) 

(BIO) 
/sinh [P(Z2-Zi)] 
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